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Abstract 

We present numerical evidence for a transition between the Klebanov-Strassler 
background and a solution describing a black hole in the class of cascading so- 
lutions in the chirally restored phase. We also present a number of properties 
of this solution, including the running of the coupling constant, the viscosity to 
entropy ratio and the drag force on a quark moving in this background. 



1 Introduction 



The AdS/CFT correspondence provides an alternative way to study the strongly cou- 
pled regime of gauge theories via gravity duals. The original statement of the AdS/CFT 
correspondence identifies J\f = 4 supersymmetric Yang-Mills with IIB string theory on 
AdS^ X [D 12] • At finite temperature the gravity side is described by nonextremal 
D3 branes and the qualitative matching of the properties was one of the key observa- 
tions in the understanding and eventual formulation of the AdS/CFT correspondence 
[3]. Other interesting aspects of finite temperature theories, as seen by the AdS/CFT 
correspondence were discussed by Witten [3j. In particular, the Hawking-Page phase 
transition in the gravity side was related to the confinement /deconfinement transition 
on the field theory side. 

In a remarkable series of papers [5], El [71 El E] Klebanov and collaborators carried 
out a program that concluded with a supergravity background that is dual to a con- 
fining field theory with chiral symmetry breaking. This background is known as the 
warped deformed conifold or the Klebanov- Strassler (KS) solution. In the context of 
the gauge/gravity correspondence it is natural to ask about the possible deconfinement 
transition for the dual field theory at finite temperature. The finite temperature phase 
of such theory is certainly very interesting and has been tackled in various papers in- 
cluding [m nil |13]. In particular, reference [13] constructed a perturbative solution 
whose regime of validity is restricted to high temperature and small value of the flux: 
Jj,^ F3 = P. Further improvements to this solution were presented in [T3|. Knowing 
the solution only asymptotically in the radial coordinate and for a specific regime of 
parameters prevents us from extracting the full thermodynamics and from being able 
to understand possible phase transitions. In particular, in order to answer questions 
that require knowing the solution for all values of the radial coordinate (for instance, 
those questions involving the action), it is important to have the full solution which 
was found numerically in |15j . 

The thermodynamic aspects of strongly coupled field theories have gained a lot of 
attention recently, largely stimulated by RHIC. Some of the thermodynamic aspects 
discussed in this paper have been considered in simpler supergravity backgrounds in 
[T9I [20I EH [22I [23]. The literature is by now very extensive, a representative sample 
of the different questions that are being pursued can be found in the pages of the 
Perimeter Institute, where part of our results were preliminarily presented [24j . 

The paper is organized as follows. We review the construction of the cascading black 
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hole in section[2l In section[3]we present our main result which is the transition between 
the KS solution and the cascading black hole. Section H] discusses various properties 
of the cascading black hole, including the running of the coupling constant at finite 
temperature, comments on the viscosity bound and the drag force. We conclude in 
section \5\ with a brief discussion of our main results and some open questions. 

2 Review of the cascading black hole 

To understand the Ansatz for the finite temperature solution we start from the con- 
formal case where the supergravity background is AdS^ x T^'^ and the field theory 
was discussed in [5]. The five dimensional manifold denoted by T^'^ is parametrized 
by coordinates (V', 6'i, 0i, 6^2, 02)- The construction of the Ansatz follows directly the 
one presented originally in [12] . For the metric we consider a generalization consistent 
with the f/(l) symmetry generated by ■i/'-rotationj^- The Ansatz in question depends 
on four functions (x, y, z, w) of the radial coordinate denoted by u: 

ds^ = e'"(-e-^"rfX2 + e^'^dXidX') + e'^'dsl , (2.1) 

where 

dsl = e^^'Mu^ + e^\dM^f , 
{dM,f = e-'-el + e^-iel+el+el + eD^e^-dsl (2.2) 
The Funfbein is: 

= i;{dip + cos9id(t>i + cos92d(j)2) , eg, = -^d9i , e^^ = sin 6'j(i0j . (2.3) 
o v6 V 6 

The qualitative meaning of the metric functions x, y, w and z was explained in [12], [15]. 

The Ansatz for the p-form fields is as in the original Klebanov-Tseytlin (KT) solu- 
tion [H]: 

Fs = Pe^ A (e^i A - A 6,^2) , 
B2 = f{u){eei A e<^i - ee^ A e^J , 

Fr^ = T + *T , T = K{u)e^ A ee, A e^, A ee^ A e^^ . (2.4) 

^In the gauge theory [5] this symmetry is identified with the U{1)r. Restoring this symmetry 
at high temperature is understood as chiral symmetry restoration [13| . Therefore, the solution we 
constructed in [TS] corresponds to the deconfined phase due to the presence of a horizon but also to the 
phase with the chiral symmetry restored due to the fact that U{1)^ is a symmetry of the background. 
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Note that the form of F3 is such that it describes a constant flux along a 3-cycle, that is, 
J^^ F3 = P. In some other notation this is called the number of fractional D3 branes. 
The Bianchi identity for the 5- form, d* = dF^ = A F^, implies 

K{u) = Q + 2Pf{u), (2.5) 

where Q is a constant. P and Q are known as M and in the standard literature (see 
review [TU]). 

Thus, in the presence of 3- form flux (P 7^ 0), the flux of F5 varies with the radius. 
The fact that K{u) depends on the coordinate u is very novel and has interesting 
physical implications. Phenomenologically, a very attractive property of this class of 
supergravity solutions is that it encodes the logarithmic running of a combination of 
gauge couplings in field theory. It does so via a varying B2 field which is compensated 
by a constant F3 flux through a 3-cycle. The five-form, which is constant in most 
solutions, varies according to the Bianchi identity dF^ = H-^ A F3 and generates a 
varying flux as J^^, F5 depends on the radial coordinate. The supergravity solution 
therefore has varying flux. This varying flux was interpreted in |9j as the dual of a 
Seiberg duality cascade, coining therefore the term of cascading solutions. 



2.1 The system of equations from reduction to 1-d 



An efficient way to derive the system of type IIB supergravity equations of motion is to 
start with the 1-d effective action for the radial evolution which follows from the 10-d 
action. This approach is convenient because to establish the transition we compare the 
actions. So, a simple action is an advantage. 

The type IIB supergravity equations of motion follow from the action 



S 



10 



12 



h'^idCf - ^e*(9C2 - CdB,f - ^F^ 



{dB^). 



12 
^diB 



,ei^CidC2d -B2 + 



-2-4!-(3!)2' 

F5 = dCi + 5{B2dC2 - C2dB2), 



supplemented with the on-shell constraint F5 = *F5. 
For the metric (12.11) 



G 



,Wy-2z 



108 



sin 61 sin 62 



(2.6) 
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J d}\VGR ^ C,.h.^ j dn[5i/'2-3x'2-2z'2-5w'2 + e^^(6e-2"'-e-i2«,)] ^ ^2.7) 

where Ci,,h. = 64 7r^/?b./i V3. The inverse temperature of the black hole is represented 
by Pb.h. and V3 is the spatial volume in the and directions. 

The matter part is 



-Cb.h.^ jdu^ $'2_^2e-*+^"-^^-^"'/'' + 2e*+'^+'^+^"'P' + e^^(Q + 2P/)2 . (2.8) 
The 1-d effective lagrangian is: 

L=T-V , 

T = 5y'^ - 3x'2 - 2z'^ - 5w'^ - ^^'^ - ^e-^+^'-^v-^"" f'^ , (2.9) 

V = -e^y{Qe-^-' - e-12-) + 1 *+4.+%+4«,p2 ^ IgS.^g ^ 2P/)2 , (2.10) 

4 8 

supplemented with the "zero-energy" constraint T + V = 0. 

This is effectively a classical mechanical system. The simplest equation is for the 
nonextremality function x (note that it breaks Loretz invariance in the {X^, X^, X^, X^) 
plane): 

x" = , X = au , a = const. (2-11) 

The reason for such a simple equation is that, as explained in [12], it does not appear 
in the effective one dimensional Lagrangian, except for its kinetic term. 

The other functions y, w, z, f and $ should satisfy a coupled system of equations: 

lOy" -Se^yiGe-^"" -6''^"") + ^" = 0, 

4z" - (Q + 2P/)2e«" - e-*+4--42'-4«'(/'2 + ^2<S>+8y+8n^ p2^ ^ 

^g_$+4^-%-4«, j/y _ p(g ^ 2P/)e«" = 0. (2.12) 
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The solutions to this system are subject to the zero-energy constraint T + V = 0, i.e. 

^1 $+4.+4y+4«,p2 ^ IgS.^g ^ 2P/)2 - Sa' = . (2.13) 

4 8 



This constraint was used in [15] as a criterion for quahty of the numerical solution and 
it was established that for various known solutions it was of the order 10^^^ where 
the reliable accuracy was shown to be 10^^° (see figiH and section 12.31 for a detailed 
description) . 

The function y amounts to a choice of the radial coordinate. This is relevant for 
understanding the dimensions of all quantities. Note that in particular, the system 
(]2.12p and the constraint ( 12.13^ are invariant under — > Lqc^ and u Lq'^u if we 
assume that Q L^Q, P ^ L and / L^f. We therefore express all 

dimensionfuU quantities in units where Lq = 1. 

Standard regular non-extremal D3-brane solution: To develop our intuition, we 
present the nonextremal D3 brane of the system fl2.12p in the radial u-coordinate. It 
takes the following form: 

^ ^ g4z ^ ^ g4x ^ g4a« _ (2.14) 

sinh4aM ' gsinh4au ' 

Note that near the horizon {u ^ oo) 

y = y,-au + ^e"®"" + 0(6"^^"") , z = - au + ^e"^"^" + ©(e"^^'^") , (2.15) 

1 1 1 , 2a , , 

y^ = -ln2a, = -In — , (2.16) 
4 4 g 

A more recognizable form of this solution is given by: 

ds^ = h~^/^{gdXl + dX.dXi) + h^%-^dp^ + p^{dM^Y] , (2.17) 

r, — ^-sx _ 1 2a 4 2a -az~ax _ 9 (o 

9-^ -^--T' P - /^-e --7- l^-^^J 



Let us clarify the relationship between the radial coordinate u and the more stan- 
dard coordinate p. Using 2a = pg^ "we have 



,„=*fi-§y'. (2,19) 



V P^ 



This can be integrated to 



Pi 



Note that in the domain of po ^ P < oo we have that u ranges in < m < oo. For 
large values of p we have that u ~ l/4p'^. The position of the horizon which is finite 
in the p coordinate becomes infinite in the u coordinates. 

2.2 Universality of regular nonextremal D3 brane horizons 

For the general solution of the system fl2.12p . in the u-coordinates the area of a surface 
defined by a horizon at m = constant is 

A = Vujr,e^^{-2z + ?,x + by). (2.21) 

Given that the equation of motion for x has the general solution x = aw we are forced 
into the following situation. If the horizon is at n ^ oo, as is the case for the non- 
extremal D3 brane solution, then for the area A to be finite we need the following 
asymptotics for z and y: 

z^aau + z^:, y—>-Pau + y^, (2.22) 

with the condition that 

3 -2a + 5/5 = 0. (2.23) 

Note that the regular nonextremal D3 brane corresponds to a = f3 = —1. The main 
claim is that: The existence of a regular horizon fixes the asymptotic behavior of the 
metric coordinates x, y and z near the horizon. Eg. (12. 231) was used in [15J as a criterion 
of quality of the numerical solution near the horizon (see table [1] for the precision of 
this criterion for the solutions we discuss here). 

Similarly, one can obtain an expression for the temperature. Namely, the relevant 
part of the metric is 

ds^ = e^'-'^dr^ + e-^'+'^^du''. (2.24) 
We introduce a new radial coordinate as: 

p = e"-^'". (2.25) 



We can now rewrite the metric as: 



ds^ = 

z' - 3x'V 



dp^ + p" (e2^-5^-3-(z'-3a) dr)' 



(2.26) 
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Note that, again, the requirement of finite temperature fixes the large-u asymptotic of 
various metric functions to he 2z — 5y — 3x —>■ constant. Requiring the absence of 
conical singularity we find that the temperature defined as the inverse of the period is 



where we used that near the horizon the asymptotic form of z and y is given by (12.221) . 
2.3 Numerics 

In Ref. P3] we gave a detailed description of the construction of solutions to the sys- 
tem fl2.12l) which contain regular black holes in the infrared (large values of u) and a 
cascading behavior in the ultraviolet (small values of u). We will summarize here our 
approach. For more details on the numerical method and the understanding of the 
numerical output, we refer the interested reader to Ref. |15j. 

To solve the equations for the metric functions we use the seventh-eight order 
continuous Runge-Kutta method. Thanks to its adaptive scheme, this method provides 
a great control upon the output accuracy. We tested and fitted the numerical procedure 
using a number of known analytical solutions, some of them with fixed values of the 
parameters (for instance, the non-extremal D3 black hole, with P = 0, and the KT 
background, with a = 0). Analyzing the sensitivity of the solutions to the variation of 
these parameters, we found that for a method tolerance of 10~^^, an error of 10^^*^ can 
be safely regarded to be negligible, what allowed us to set our numerical 'zero' to this 
last value. 

Using the non-extremal D3 solution (P = 0), we set the boundary conditions for our 
ten variables at the value of u corresponding to the 90% of the distance to the horizon. 
A correction is introduced to enforce the solution with P 7^ to satisfy constraint 
(12.131) in this boundary. 

Next we integrate backward and, for a given value of a and Q, we find the value of 
P such that the space is complete, i.e., Usmg ~ (see section [2^ and table [1]). We use 



T 



a — 3 
27r 



a e 



.2z,-5y, 



(2.27) 



fi + Pe 




(2.28) 



to check that this singularity is of the type given by 




(2.29) 
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a 


P 


""90% 


3 - 2a + 5/? 






1 


4.26906025 


2.19361453 


-9 X 10-15 


3 X 10-1° 


0.00001257 


0.5 


2.87343889 


4.387229 


-7 X 10-15 


3 X 10-11 


0.0000065 


0.1 


1.094389139 


21.9361453 


-2 X 10-15 


4 X 10-1° 


0.00000002 


0.01 


0.316430028 


219.361453 


10-15 


2 X 10-9 


0.0000000002 



Table 1: Values of the parameters used for the four cases presented in this manuscript. 
Q = Lp = 1 was used everywhere. 
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Figure 1: The Hamiltonian constraint for the cases in tabled! 

and that the remaining fields are analytical at Using- In this criterion, the subindices 
indicates that these are coefficients in Taylor series for the corresponding field at Using- 

We then integrate the system forward and identify numerically the presence of a 
horizon. One of the criteria used is that given by Eq. fl2.23p . 

In table [1] we present the values of the parameters and the results for the infrared 
and ultraviolet criteria of quality for four cases that are going to be used as examples 
throughout this manuscript. 

We also recall that system (12.121) was modified in such a way that the solutions must 
automatically satisfy constraint fl2.13p . This allows to use the constraint as a dynamical 
criterion of the quality of the numerical output for any value of the parameters. Indeed, 
all our solutions are warranted to satisfy constraint fl2.13p within the accuracy given 
by our numerical 'zero'. This can be verified in figiU For all the four cases in the table, 
the T + V = constraint is satisfied in the whole range of the radial coordinate. 
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2.4 UV asymptotics and completeness of the solution 

Let us first discuss the completeness of our solution. The statement that completeness 
of the space requires Using = was made explicitly in (see section 4.2 there), 
however, the details were not shown. Here we provide those details since they are 
crucial in clarifying the number of parameters of the solutioij^. 

The asymptotic form of the relevant part of the metric, as obtained in [15], is given 
for z{u) by Eq.f l2.29i) and for y{u) by exp{Ay) = l/Au. The first expression was widely 
discussed in [15] and we refer the reader to that paper for details, including the fittings. 
The expression of y is ubiquitous in the region of small values of m. It is valid in the 
asymptotic region for all the solutions we know: KW [5] , nonextremal D3 branes [12] , 
KT [8] and KS [9]. The best fit analysis presented in [15] indicates that this seems to 
be also the case for our numerical solution. 

To discuss completeness we consider an outgoing null geodesic toward the asymp- 
totic regions of small u. The corresponding effective Lagrangian is: 

£ = _e2-6x^2^e-^^+i°V. (2.30) 

With this information the affine time for a geodesic to reach Usmg satisfies: 

d\ = ^u-^/^du. (2.31) 
E 

Note that this is independent of the precise form of the z function in the metric. The 
structure of the function z enters only through the integration limits. Namely, the 
above equation can be integrated to yield 

,2.32) 



\^sing ^i,. , 

Here Um should be viewed as a cutoff, beyond which we do not know the precise 
formula for y. However, we assume that this geodesic originates somewhere near the 
black hole horizon. Note that the value Using enters just through the integration limits, 
since the full solution does not exist for u < Using- As can be seen from (12.321) the 
affine parameter is infinite only for Using = 0. We need an infinite affine parameter 
to guarantee completeness of the solution. Thus, the solution is complete only for 

"^sing . 



^This point has been an important source of discussion with A. Buchel and we hope that this 
section clarifies that Using = fixes one of the three seemingly free parameters {Q, P, a). 
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The number of relevant parameters 

Thus, superficially, our solutions depend on three parameters, a, P and Q. However, 
the condition Ugmgi,- ■ ■ ',cl,P,Q) = imposes a new constraint on the solutions that 
effectively reduces the number of parameters to two. 

For instance, fixing Q = 1, for the values of a given in tabled! it is found that the 
expression, 

Us^ng = fe"^' + 60 + hP + 62^' + " " " , (2.33) 

fits very well the numerical output. Note that the first term is clearly related to the 
strong energy scale as indicated by the KT formula, that is, it is related to the radial 
position where the KT warp factor vanishes; of course, there are many corrections as 
the KT warp factor is not the full answer. To illustrate the origins and robustness of 
the expression (12.331) . let us first discuss a specific point in the graph [2] to highlight the 
importance of the first term in the above expansion. Namely, for a = 0.5, a best fit to 
the value of Ugmg without the exponential term (i.e., setting b=0) yields the following 
result for the coefficients bf. 

{0.05854539, -0.25108172, 0.14454983, -0.03078276, 
0.00347435, -0.00020362, 0.488 x 10"^} , 

with a goodness of fit on the order of 10~^. It is important to stress that for each term 
we verify that 6j < fej+i, which indicates convergence. Now, introducing the exponential 
term, the result is 6 = 5.49169633, 

{-0.09506874, -0.09298412, 0.07882655, -0.01663906, 
0.00180981, -0.00010186, 0.234 x 10"^} , 

with a goodness of fit on the order of 10~^. The ratio is also much better 

behaved. In general the exponential term seems to be crucial for the form of Ugmg- 
Thus, for (a, Q) fixed, by setting Using = in expression fl2.33p . we are this way 
selecting a given value of P. A similar analysis was carried out for the other values of 
a in table [1] and another ten cases not included in the table but shown as crosses in 
figl21 The curve in this figure represents the pairs (a, P) that guarantee Using to vanish 
for Q = 1. It was obtained as the best fit solution, 

P = 4.26846923 a°-™^^^ + 0.02429978 



10 



p 



2 



4 



3 



1 




0.2 



0.4 



0.6 



0.8 



1 



a 



Figure 2: Values of P and a that guarantee Usmg to vanish for Q = 1. 

with a goodness of fit of 0.00124, which is satisfactory if the small amount of data used 
is taken into account. 

Summarizing, for any values of two of the parameters determining our solution, 
the value of the third one is fixed by imposing completeness of the solution in the 
asymptotically UV region. 

3 The Hawking-Page transition 

There is mounting evidence for the transition we are going to describe here. In par- 
ticular, similar transitions have been established in some simpler holographic models. 



In our case, we assume that the deconfined phase is dual to the cascading black 
hole reviewed in the previous section, while the confined phase is dual to the warped 
deformed conifold, i.e., the Klebanov-Strassler background [9]. The relevant expressions 
of this well known solution were compiled in the appendix in terms of the coordinates 
and notation we use in this paper. 

To accurately compute the actions we used an adaptive 3-5 Simpson's quadrature. 
This implies 2"~^8'^ evaluations of the Lagrangian, where n is the number of iterations 



including [29l EQl EH |32] . 
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needed to reach the desired precision (set here to 10^'^) for the integral. For the 
cascading black hole the Lagrangian is given in terms of the numerical solutions of 
system f l2.12p . We gave priority to the output accuracy over the computational effort. 
This means being able to obtain the actual numerical solution of (12.121) at every point 
required by the quadrature, instead of finding it by interpolating between previously 
calculated solutions at a given set of points. The whole computation of the cascading 
black hole action usually takes more than 12 hours in an ordinary desktop computer. 
Compared with that, estimating the action of the KS background is a quick task, since 
the evaluation of the corresponding Lagrangian only involves the numerical integration 
of (1A.14I) which for large r can be accurately estimated using expression flA.lTD . Taking 
the above into account, we decided to compute the cascading black hole action for the 
four sets of parameters given in table [Hand compare the results with the output of the 
KS action as we vary e. 

So, we calculated the cascading black hole action as given by Eqs. (l2.7p and (12.81) . 
without considering the numerical factors common to both actions. The integration 
was done between = 10~^ and ugQ%. Following the order in table [Tj the four 
results are {738451, 762361, 733711, 804269}, and their logarithms are shown in figlHas 
horizontal lines. To leading order the asymptotic values of the action are independent 
of a. This can be seen by directly evaluating the pieces of the action as given in (12. 7p 
and (12.81) . For example, the leading contribution from (12.71) comes from z'^ and after 
the integration is proportional to 1/ur. Similarly, from the matter part (12.81) . the 
leading contribution comes from the last term which is independent of a; this term also 
contributes a term proportional to 1/ur. Note that in terms of the standard radial 
coordinates, this corresponds to i?^ where i? is a UV cutoff radius. This is the typical 
volume divergent behavior of AdS^ as shown in [U [211 ED] ■ 

To compare these results with those for the KS action is necessary for both back- 
grounds to have the same physical temperature. This is implemented by requiring that 
the physical perimeter of the temporal directions match at the point of comparison: 

f^^^^z(ua)-Mua) = (3j,sh-^l\TR). (3.1) 

This equation allows to determine I3ksi since for the black hole, /^^./i. is determined 
by the absence of conical singularities (see equation (I2.27P ). However, we still need a 
way to relate tr and ur. As can be seen from various examples discussed in [121 [15] . 
in the UV the coordinate u is related to the standard radial coordinate as m ~ l/4r"^; 
this is also true for the nonextremal D3 brane (see equation (I2.20p ) as well as for the 
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KT solution, the b-deformed conifold ^2J, KS and the resolved conifold [18J. The 
relationship between the radial coordinate r and the standard conifold radius is 

e^/^ sinhr 
= 21/63V2 (sinh2r-2r)i/3^^- ^^'^^ 

In the limit of large values of r we obtain: 

- - (3-3) 
Thus, in equation (13.11) we use that 

tr^^ \n{e-^u^^) - 0.954771252 + 6tr. (3.4) 

The correction 6tr is required to account for the limitations of the numerical analysis 
to deal with actual asymptotic values. It is a function of a and e and we have found it 
to have the approximate form, 

Str = C(a)ei/^ . (3.5) 

Our strategy to establishing the existence of the transition is as follows. We assume 
that the transition takes place at a value of e such that a~^^^e = 1, then we find the 
C{a) that makes both actions equal. Following the order in tabled! we obtained that 
C{a) = {-2.27, -3.264, -5.85, -9.8095}. We note that neither the factor C{a) nor 
the power 1/5 are solutions of a best fit problem but just values that give satisfactory 
results for all cases analyzed here. Our key observation then reduces to establish that 
the difference of actions changes sign appropriately around this point. 

Let us address a technical issue. A problem to deal with here is the singularity 
in the KS background for Q ^ (see in the appendix subsection lA.l.ll for details). 
Fortunately, although conceptually an important point, the difference between the 
KS actions with and without Q are qualitatively distinguishable only in the infrared 
regime (see for instance figures [3]). In the UV regime, the leading asymptotic behavior 
is governed by P and not Q. 

The IR behavior is easy to understand. In this region, almost at the same value of 
r, the Lagrangian with Q = has a zero and the corresponding Lagrangian with Q ^ 
has a minimum. We therefore switch at this point from the integrand with Q ^ to 
the one with Q = as a regularization of the procedure. From the quantitative point 
of view, for Q = the infrared regime gives a negative contribution to the integral 



13 




Figure 3: The KS Lagrangian in the IR and the UV for e = 0.001 and e = 1 
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Figure 4: The Hawking-Page transition. The horizontal hnes represent the values of 
the cascading black hole action for the four cases analyzed in this paper. The curves 
show the e-dependence of the corresponding KS actions. 

which increases slowly with e. Thus, for a fixed a, if for small e the cascading black 
hole action dominates the difference of the actions, this is a reliable result because the 
negative contribution to the integral is negligible. On the other hand, if for large e the 
KS action dominates, that is certainly true, because it does so in spite of the negative 
contribution to the integral. 

Finally, the results of the actions comparison are presented in figJH Curves represent 
the variation with e of the KS action. Each curve with a given color must be compared 
with the horizontal line of the same color, i.e, with the cascading black hole with the 
corresponding temperature. As we can see, in each of these cases, the deconfined phase 
dominates for small values of e, while the confined phase dominates for large values 
of this parameter. The larger the value of a, the larger the value of e for the phase 
transition to take place. 

A limitation of the approach described here is that it works for relatively small 
values of e, implying relatively small values of a. Otherwise, the correction given by 
(13. 5p might stop to being small. Moreover, the difference between the KS solutions 
with (5 = and Q 7^ could contribute as well. Our calculations indicate that the 
transition can be detected for values as large as e = 10 (a ~ 1). 
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4 Some properties of a confining field theory at finite temper- 
ature via gauge/gravity duality 



4.1 The running of the gauge couplings at finite temperature 

One of the most interesting properties of the cascading solutions is that they encode 
the running of a gauge coupling. In particular the matching of the supergravity modes 
with the two gauge couplings takes the form [TU] : 
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47r2 
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9l 


9l 


9se^ 


2-na' 




TT 



(4.1) 



For the KT and KS solutions one finds a result that exactly matches the field theory 
value of the running coupling for a gauge theory with SU {N + M) x SU{N) gauge 
group. 

For the cascading black hole solutions we can assume that a similar correspondence 
holds and evaluate B2. The integral of B2 over the appropriate two-cycle is proportional 
to the function f{u) in (12. 4p . Since our solution is numerical we need first to find out 
what is the behavior of f{u) in the UV regime. If we assume that as u — > 0, 

/(n)^-|pe*(")ln(aw) + /,, (4.2) 

where fr is a constant, then criterion (I2.28P fixes 



B = lim 4 



.0 $1 + i . 

^ u sing 

The pole in this expression is expected to be cancelled by the logarithmic singularity in 
Eg. (14.21) . Substituting B in (14.21) . we compare the behavior of the resulting expression 
with the numerical solution for f{u) near Using- Denoting their difference as A/, the 
obtained results are presented in figure ([5]) for all four cases in tabled! 

It can be seen, in the UV regime, that up to corrections well below 10^^*^, fr is 
indeed a constant. The approximated values of corresponding to the analyzed cases 
are {3.22186396, 2.22309019, 0.99140974, 0.41584279}. 

Thus, the running of the difference of the inverse squares of the coupling constants 
have several interesting properties. First, it is proportional to P which is the parameter 
breaking the conformal invariance. More importantly, we have the following 
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u 




Figure 5: Behavior of A/ for a = 1, 0.5, 0.1, 0.01. Since it is clearly below our 'zero' 
of 10~^°, this shows that equation (14. 2 p represents very accurately the form of f{u) in 
the UV-limit. 
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Hau) ~ In ~ In . (4.3) 

In the above expression we have used that 2a = pg, where po is the position of the 
horizon in the standard nonextremal D3-brane solution. Further po = TttR"^ for nonex- 
tremal D3 branes and finally u ^ 1/ Ap^ in the uv region. 

Thus, we conclude that the gravity dual of the cascading black hole has 



It is worth stressing that the high temperature behavior of the gauge coupling in 
Yang-Mills theories depends on the renormalization schemes at finite T. There also 
seem to be some discrepancies between the imaginary and real time formalism. In a 
sense our result should be interpreted as the string theory prediction and it coincides 
with the general perturbative reasoning [23 EH]- 

4.2 The viscosity bound 

In this section we simply note that the solution [IS] satisfies the condition of [27] which 
implies the viscosity bound. Namely, we find that for the metric (12.11) the Ricci tensor 
satisfies ^ 

R%-R\=Ae^^-'''y (4.5) 

This expression is identically zero due to the equation of motion of x. 

In principle, it has been argued that the above expression is sufficient to guarantee 
the viscosity bound [271 [28]. We will nevertheless proceed to check the bound explicitly. 
The main reason being that our background is cascading and that might involve sub- 
tleties not considered in previous analysis and also that [27j used some nonconvariant 
counterterms. Since our solutions satisfy the criterion (12.231) . we obtain for the shear 
diffusion constant, 

D 



la — 31 1 



167r T 

Assuming the black hole temperature to be given by (12.271) . for the viscosity-entropy 
ratio, 

rj |a — 3| 
s IGvr 

In Fig we can observe that, for all the four cases in table [H a — — 1 from below 
as u approach the horizon. In this figure the values of u were conveniently rescaled for 
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Figure 6: Behavior of the parameter a as m — Uhor- 



each case using the ratio between the different values of a. The result seems to indicate 
that indeed the inequality is saturated. We must note here that, in all the cases we 
analyzed, the value to which a converges is actually slightly bigger than —1. For these 
four cases we obtained {-0.999999937, -0.999999946, -0.999999977, -0.999999997}. 
Nevertheless, since the difference with respect to —1 is too close to our numerical 
'zero' we regard this effect as a numerical artifact due to the proximity to the horizon. 

4.3 Drag force 

In this section we apply the analysis of [Ml [35] to discuss the drag force on a quark 
moving in a strongly coupled plasma. We emphasize the general conditions for a su- 
pergravity background to admit the motion of a classical string that can be interpreted 
as a drag force on a moving quark. Our analysis mirrors the generality discussed 
in [36], however, here we do not assume AdS asymptotics, in particular we aim at 
understanding cascading backgrounds. 

We consider here a five-dimensional subspace of a supergravity background dual to 
a field theory. 



where the metric components Gmn are functions of the radial coordinate u only. Fol- 



ds^ — — GqqXq -\- Gxxi^dX'^ -\- ~l~ dX^^ -\- G^udu^ ^ 



(4.6) 
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lowing [3ll [35] , we assume the world-sheet to be embedded as t = Xq and a = u and 
we allow for 

X, = vt + au), (4.7) 

which means that the end of the string is moving with velocity v. In the field theory side 
we interpret the end of the string as a quark moving with velocity v. The Nambu-Goto 
action is 

S = ^j dtda^[GooiGuu + - G,,GuuV^]. (4.8) 

Since the action does not depend explicitly on a, one has that the conjugate momenta 
1^ is a constant, where L is the Lagrangian, 

= _ = — "^"^^ (4.9) 



.1 n„„(n J- n f^i^ 



Goo{Guu + Gxx^' ) ~ GuuGxxV 



2 



It can be rearranged to obtain 



A key simplifying observation is that, since ^' cannot be imaginary, we need both 

expressions in numerator and denominator to flip signs simultaneously. This fixes 11^ 
to be 

n| = G^o{x)GxxM. (4.11) 
Here, denotes the radial coordinate satisfying 

GooM = Gxx{u,y. (4.12) 

The rate of change of momentum is calculated to be 



dp 

where 



^ = V^T\,, (4.13) 



= -^Gx,ug'"'<PaX''. (4.14) 

Here, Gij denotes the metric in fl4.6p . while gij denotes the induced metric on the 
world-sheet. After some algebraic simplifications, we obtain 

^ = -^n,. (4.15) 

dt 2TTa' * ^ ^ 



20 



We can rewrite of above equation as 

^ = -kPi. (4.16) 

if condition (14. 12p dictates 11^ to be proportional to pi. In the above equation, k is the 
drag force constant. 

Conditions for drag force 

Let us discuss the conditions under which we can observe a drag force proportional to 
the momentum. We need to solve for a radial position m^, such that 

= Goo(m*)/G..(m*). (4.17) 

The expression for the canonical momentum is 

= vGUuM)- (4.18) 
So, the question for a drag force essentially boils down to verifying that 

GUu.{v)) = -7=L=, (4.19) 
V 1 — v'^ 

where 6 is a constant. 

Two comments are in order about the relevant energy scales and the speed of the 
quarks: 

• The nonrelativistic motion of the probe corresponding to small v is localized near 
the points where Goo vanishes. This corresponds, generically, to the horizon. In terms 
of energy scales, it corresponds to the infrared region. 

• Relativistic velocities v close to the speed of light correspond to radial positions for 
which Goo ~ Gxx- Generically, this is the asymptotic region, which in terms of the 
field theory is the ultraviolet region and corresponds to the near conformal limit. 

Non- extremal D3 brane 

In the case of non-extremal D3 branes we have that 

Goo = u'(^l-^y = u\ (4.20) 

The equation that determines the position Vi^ clS 3i function of the velocity of the end 
of the string is thus: 

2 Gqq ( ^ ui 
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Thus, we find 



We can now easily verify tliat 



ul = ^=^. (4.22) 



ul 



GUu.{v)) = ul = ^=^ (4.23) 
V 1 — V 

Tliis expression coincides witli tlie criterion ( 14.19p . Tliis case was discussed in Refs. 



Cascading black holes 

Here we simply adjust the computation of the drag force presented, for example, in 
[211 [35]. The calculations presented here were, in particular, performed in [371 EH] 
without using the full metric of the cascading black hole. The AdS part of the metric 
in consideration in string frame is 

ds^ = e^l\^\-e-^''dXl + e^'^dXidX') + e^^'y-''^ du^\ (4.24) 

We consider the worldsheet to be along t = Xq and a = u directions. Following the 
general arguments at the beginning of this subsection 

= v^e*+2-(u(a;(i;))) (4.25) 

Here, u is chosen to satisfy v = e~^^*^"^ for a given v. This relation is to be viewed as 
an equation for the radial coordinate u, it defines a particular value Up. We need to 
evaluate relation fl4.25p at u = Up and, in turn, find H^ as a function of v. Further, to 
establish the existence of a drag force we need to verify the existence of a 6 such that 

He = b^=L==, (4.26) 



y2 



which thus implies that H^ = and directly leads to a drag force parameter using 
equations (14.151) and (14.161) . With this aim, equating equations (14.251) and (14.261) we 
obtain, 



b{u) = v^e^+^Vsinh (4 aw) 

v{u) = e"'*^". (4.27) 
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Figure 7: Behavior of the drag force parameter 6 as a function of the quark velocity v. 
First, note that in the IR regime {u oo), 



where we used the asymptotical expression fl2.22p for z{u). In the pure D3 case, where 
= and a = —1, we have that b = exp{2z^). Recalling that, as shown in 
subsection 14.21 for our solution a ~ —1, and noting that the dilaton is very weak 
near the horizon (for instance, for our four cases we obtained ^ {—3 x 10^^, —2 x 
10^^, —1.2 X 10^^, —1.4 X 10^^}), we see that for the cascading black hole b ~ constant 
in the IR-limit. 

On the other hand, knowing that as m ^ 0, z{u) behaves as described by Eg. (12.291) . 
we obtain. 



implying that b is finite in the UV-limit. 

This way, we have shown analytically that when the string explores any of the two 
asymptotic regions, we verify the existence of a drag force. For a given temperature the 
variability of b from one regime to the other is mainly related to the dilaton variability. 
In figl7]we show the drag force coefficient b as given by Eqs. (14.271) for all values of v, 
that is, for all velocities of the external quarks. The curves represent the four cases 
in table [H It is worth stressing that in the relativistic regime there is a clear velocity 
dependence in the coefficient b. However, for the range of values of temperatures that 



b 



(4.28) 




(4.29) 
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we have explored it, it does not seem to be universal. In particular, the fig 17] shows 
that in the regime of relativistic velocities some curves have infiections upwards while 
other bend downward. Clearly, it seems to depend on the temperature and deserves 
further study. 

5 Discussion and outlook 

In this paper we have established numerically the existence of a Hawking-Page phase 
transition between the supergravity solution of [15] which represents a black hole in a 
cascading background and the Klebanov-Strassler solution. 

We also studied some properties of the deconfined phase given by our cascading 
black hole solution. For instance, we have found that the difference of the inverse 
squares of the coupling constants is proportional to P (which is the parameter breaking 
the conformal invariance) and that it runs as ln(A/T). 

We also verify that our solution seems to satisfy the well-known bound for the 
viscosity to entropy ratio. 

Finally, we have explicitly shown that the universal form of the horizon in this 
theory guarantees the existence of a drag force for a fundamental string moving in this 
background and in the dual field theory that can be interpreted as jet quenching. The 
near horizon region describes nonrelativistic motion of the quarks and since the near 
horizon geometry has been shown to be universal we conclude that the nonrelativistic 
drag force on the quarks is also universal in the class of cascading theories. In the 
relativistic region we also found a well-defined drag force parameter. Numerically, we 
observed a velocity dependence for this parameter in the relativistic region. 

It is important to stress that most of the calculations performed in this paper 
requires knowing the solution for all values of the radius and that an analysis solely 
based on the asymptotic near the horizon is not reliable. 

Let us conclude with listing a number of possible directions. Cascading theories, 
along the lines of the particular model described here are generic in the AdS/CFT. 
In fact, we believe that there might be infinite classes of them based on some results 
presented in [5^ HO] . The fact that the cone over Y^'*^ does not admit complex deforma- 
tions seems like a deterrent for a construction directly mimicking the warped deformed 
conifold of Klebanov and Strassler but a more general solution where the fiuxes play 
a crucial role are not ruled out. Turning on temperature of these solutions should not 
change substantially the asymptotic UV form obtained in 1^ . 
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It would be interesting to better understand the charges in the comparison between 
the cascading black hole with the KS background. Basically the KS solution has 
Q = 0, that is, the number of D3 branes at the tip deformed conifold equal to zero. 
As shown in [H] and reviewed in section lA.l.lj we see that nonzero Q results in a 
singularity at the origin. Another interesting direction would be the construction of 
black holes in the backgrounds with back reacted flavor such like the one described in 
Ref. pTj. In particular, some of these backgrounds should be natural generalizations 
of the construction of [15]. 

While this manuscript was in the final stage of its preparation, a new paper was 
posted where the authors claim to have found evidence for a transition in a solution 
with a black hole in the IR and KT in the UV [12]. We plan an in-depth discussion of 
the relation between our and their solutions. Preliminarily, let us note here that there 
are some key differences between our solutions. First, as far as we understand, to find 
their solution they need to numerically look for a zero of a function in ten variables, 
i.e., the norm of their 'mismatch vector'. Since this is a positive semi-definite function, 
assuming that this zero exists, this problem can be translated into finding the global 
minimum of a lOD-surface. Nevertheless, this is a case of constrained optimization, 
since the Hamiltonian constraint is also valid in ten dimensions, where it shows as 
a reparametrization of the radial coordinate. However, the authors of [42J made no 
mention to whether their numerical solution satisfies this constraint. As opposed to 
that, our solutions are forced to automatically satisfy the zero-energy constraint and 
that this certainly happen is one of the main criteria of quality we use to certify 
our numerical output. This was first explained in Ref. [15] and emphasized in this 
manuscript. Secondly, their solution seems to have exactly KT in the boundary (we 
note here that the KT solution does not satisfy the Hamiltonian constraint for a ^ 0) 
and to erase the large u singularity typical of KT (we called this singular point as 
Using), pasting the black hole metric in the IR solution somewhere in the 'middle' of 
the radial domain. Our solution extents from Using to infinity and, as explained in 
Section 12. 4[ we set Using = to guarantee the completeness of our background. 
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A The Klebanov-Strassler background 



The KS metric is of the form: 



ds'' 
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h-^''^{T)dXndXn + h^''^{T)dsl , 



dsl = '-e'/'K{r) 



the one-forms are defined as 



where 
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V2 ' 



V2 ' 



2 - 
4 e' + e^ 



(A.l) 



(A.2) 



■ sin ^1^01 , 



d9, , 



e = cos ip sin ^2*^02 — sin i/jd92 , 
= sin ip sin 92d(f)2 + cos ■ijjd92 , 

6^ = ^-0-1- cos 9id(l)i + cos ^2C?02 ■ 



The function in the metric is 



/C(r) 



(sinh(2r) - 2r) 



- 9^^1/3 



(A.3) 



(A.4) 



2^/3 sinhr 

For asymptotically large values of r one can introduce a natural radial coordinate 
= -^e'^l'^e^'^l'^ such that the metric on the conifold becomes ds\ ^ dr"^ + r^cis^i,i. 
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The matter fields are as follow: 

Ma' 



{/ A A / + d[F{r){g' A / + / A g% 



2 

^{/A/A^^(1-F)+/A^^A/F 
+F'dT^{g'^g^ + g''^g^)] , 



with F(0) = and F(oo) = 1/2, and 

S2 = ^[/(ryA/ + fc(r)5'A/], 



rfrAaVA^^ + fcVA^^) 



+^(^-/)/a(/a/ + /a^^) 

The self-dual 5- form field strength may be decomposed as F5 = !F^-\--k!F^. 
■F5 ^ A F3 ^ ^'^Y^\ {T)g' Ag'Ag'Ag'Ag^ 

where 

i = f{l-F) + kF , 

and 

★ JF5 = 4g,M^(a'fe-^^^dx^ A dx^ A dx^ A dx^ A dr ^^^^ „^ ^ . 



a = A{gsMa'fe-^/^ . 
The solution for the functions defining the matter content is: 

sinh r — r 

" 2 sinh r ' 
Tcothr - 1 
= 2 sinh r (^^^^"-1)' 

= ^T^^^I- (coshr + 1), 



2 sinh T 

T coth r — 1 
4 sinh^ T 



£(r) = /(I - F) + A;F = . (sinh 2t - 2t) . 
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The warp factor with the condition that it vanishes at infinity is 

h{T) = J^I{t) = ((7,M«')'2'/'£-'/'/(r) , (A.13) 

where 

I{t) = /"rfx^^^^^^|^(sinh(2x)-2x)^/^ (A.14) 
Jr sinh X 

Two important hmits are: 

/(r ^0) ^ao + 0(r') ; (A.15) 

/(r ^ oo) ^ 3 ■ 2-1/=^ (r~\) e'^^/^ ^ (^.16) 



4^ 

where oq ~ 0.71805. At large r the integrand becomes 



,^3.2-.i-^,„(^;_Z_j, (A,17) 

where we have used that for large radius = -^e^l'^e^'^l'^: 
A.l The action for the KS model 

A convenient way to look at this solution is through the prism of a one-dimensional 
system as discussed in [l8j. Motivated by the form of the deformed conifold metric we 
make the following ansatz for the metric 

ds^ = e^P-^e^'^dx^dx^ + du^) + [e-^^-''gl + e^+^gl + gl) + e^-^(^2 ^ ^2)j _ ^^.18) 

The ansatz for the p-forms is: 

Hs = duA [f'iu)gi Ag2 + k'{u)g^ A g^] , (A.19) 

Fs = F{u)gi A (72 A (75 + [2P - F(m)](73 A (74 A (75 + F\u)du A ((71 A (73 + A (74) , (A.20) 

F5 = ^5 + H , -^5 = A (72 A (73 A (74 A (75 , (A.21) 

K{u) = Q + k{u)F{u) + f{u) [2P - F{u)\ , (A.22) 

where F, /, k are functions to be determined and P and Q are constants. We explicitly 
ensure that the Bianchi identities for the p-forms are satisfied automatically. The 1-d 



28 



action reproducing the resulting equations of motion restricted to the above ansatz has 
the following general structure 



S = c du e 



3A'2 - \gM^"^^" - V{^) 



(A.23) 



where c = — 4|^. It should be supplemented with the "zero-energy" constraint 



3A'2-lG,,(</p)v.V' + n</') = 0. 



(A.24) 



The action (1A.23P is thus a classical mechanical action for the fields A and = 
(x, $, /, fc, F). The corresponding kinetic and potential terms in flA.231) are found 
to be: 
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(A.25) 



V{^) = ie-^P-^^ - e^f-^^ cosh y + ^e^^ sinh^ y 



ae 



(A.26) 



where K is the combination of the independent functions /, k, F given in (]A.22|) . 
The first order equations for the independent functions A, x, y,p, f, k, F, $ are: 



X 



y' = e^^ sinh y 



1 

— ( 

6 



1 

— ( 

6 



A' 



3- — y 3- Q 

1 



(A.27) 
(A.28) 
(A.29) 



f = e^+^P+^y{2P-F) , fc' = e*+^P-2^F , F' = ^-e'^^+^^if -k) , $' = . (A.30) 

The functions that we introduced in this subsection are explicitly given in terms of 
the solution of the previous subsection as: 



1 



e^/2/Csinh(r)/ii 



/2 



/3 



tanh( 



We will also use explicit forms for the derivative with respect to u 
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6 



1 

P' 



-ie^^ie^ + e"^) - e-'^'e"'" + Ke^fe"'"), 
6 



y' = e^Psinh(ln(eJ')), = -e'^Pe-^^ - ye^^'e-2^ $ = $' = 0, 
k' = Fe*e^Pe-2^ /' = (2P - F)e*e^fe2^ F' = -LzAq-^q^p, 

The expression for the action we need to evaluate is finally of the form: 

= -^e2-e-^P(3A'2 - - V), Sks = j drC. (A.31) 

A. 1.1 The KS background has no regular D3 branes at the apex 

As noted in |9j, the KS-like Ansatz allows for general solutions having K{u) with 
arbitrary Q in ( ]A.22I) . In the more standard literature Q is identified as A^. These 
solutions are, however, singular. Namely, near the apex of the deformed conifold (small 
values of r) the warp factor takes the form: 

h^^. (A.32) 

T 

An alternative way of understanding this singularity is by realizing that it corresponds 
to the freedom to add a homogeneous solution to the warp factor. Namely, a solution 
of the form 

(Vg'eg^^drh) = 0, (A.33) 



/96 

or 



dr 

K? sinh^ r 

This last statement allows to interpret the elimination of the singularity as a state- 
ment about regular D3 branes. It is possible that, as in the interpretation for the 
resolved conifold of [18J, this divergence signals a smearing of D3 branes along a three- 
dimensional subspace and therefore could also be cure, not just by turning off the 
charge Q, but also by localizing the branes appropriately as it was done recently in 
Ref.[33]. 

This situation presents a conceptual problem. It is natural to compare a field theory 
with parameters (A^, M, e ~ ^strong) with a theory with (A^, M,a Temperature). In 
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realistic cases, the scale Astrong is dynamically generated but supergravity methods are 
still far from achieving that. However, for KS we have = and for the cascading 
black hole we have Usmg = 0. The comparison is really not completely clear. We deal 
with this point while analyzing the transition in the next section. 
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